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coordination to CO molecules 
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We show that the magnetic state of individual manganese phthalocyanine (MnPc) molecules on 
a Bi(110) surface is modified when the Mn2+ center coordinates to CO molecules adsorbed on 
top. Using scanning tunneling spectroscopy we identified this change in magnetic properties from 
the broadening of a Kondo-related zero-bias anomaly when the CO-MnPc complex is formed. The 
original magnetic state can be recovered by selective desorption of individual CO molecules. First 
principle calculations show that the CO molecule reduces the spin of the adsorbed MnPc from S = 1 
to S = 1/2 and strongly modifies the respective screening channels, driving a transition from an 
underscreened Kondo state to a state of mix- valence. 



Control over the magnetic moment of a molecule and 
its interaction with a substrate is a key issue in the 
emerging field of molecular spintronics pQ. In metal- 
phthalocyanines and metal-porphyrins the metal center 
is usually coordinatively unsaturated and presents a local 
reactive site, which opens a unique possibility of control- 
ling the magnetic moment in-situ by external chemical 
stimuli [2--5J. The axial coordination of small molecules 
like CO, NO or O2 to those complexes substantially alters 
their electronic properties, which has led to a success- 
ful implementation of phthalocyanines and porphyrins in 
gas sensors [9] [10]. Studies specifically addressing the 
magnetic state are however scarce. Only recently it has 
been shown that the attachment of NO molecule to a 
cobalt-tetraphenylporphyrin (CoTPP) quenches its spin 
due to the oxidation process [6]. The general picture is 
however more complicated, as the chemical bond to the 
reactant molecule causes the redistribution of charge in 
the d-orbitals of the metal center and modifies the lig- 
and field of the metal ion. This has critical consequences 
for the magnetic ground state of the complex [7J E| . On 
metal surfaces, the formation of a new ligand bond may 
additionally alter the hybridization of molecular and sub- 
strate states, thus affecting the electronic and magnetic 
coupling of the metal ion to the substrate [3 -5 . Under- 
standing the response of these effects to the change in 
chemical coordination is crucial to gain the full control 
over the functionality of the magnetic system. 

Here, we show that the magnetic moment of a man- 
ganese phthalocyanine (MnPc) molecule (see Fig.[TJa)) 
on a Bi(110) surface is reduced when coordinated to a 
CO molecule. Using a combination of a low tempera- 
ture scanning tunneling microscopy (STM) and density 
functional theory (DFT) we find, first, that the spin of 
the MnPc molecule is reduced from 5 = 3/2 to 5 = 1 
upon adsorption and, second, that CO further reduces 
the spin of the MnPc from 5 = 1 to 5 = 1/2. The 
change in the magnetic ground state upon CO attach- 



ment is detected by the broadening of a characteristic 
zero-bias anomaly (ZBA). We interpret this broadening 
as a transition from a Kondo regime, for the bare MnPc 
on Bi(110), to a regime where charge fluctuations drive 
the CO-MnPc complex into a mix- valence state. We fur- 
ther show that the original magnetic state of MnPc is re- 
covered after tip-induced desorption of the CO molecule. 

Our experiments have been performed in a custom- 
made scanning tunneling microscope (STM) working in 
ultrahigh vacuum and at low temperature (5K). As a 
semimetal, Bismuth has a low density of states close to 
the Fermi level [TT], favoring that molecular adsorbates 
lie weakly affected by the surface [12]. We chose the 
(110) surface because it additionally presents dangling 
bonds that can eventually anchor the molecular species 
to a fixed site, creating a full commensurate layer. We 
exposed an atomically clean Bi(110) single crystal surface 
at room temperature to a flux of MnPc molecules (Sigma- 
Aldrich) thermally sublimed in vacuum from a crucible, 
and posteriorly cool it down for STM inspection. 

On this surface MnPc molecules self-assemble in 
densely packed islands (FigQa)). MnPc molecules ap- 
pear as a clover-like protrusion with a bright center corre- 
sponding to the metal ion. The adsorption configuration 
could be established with precision by resolving simul- 
taneously intramolecular structure of MnPc and atomic 
structure of the underlying substrate (see e.g. Fig.JTJc)). 
As shown schematically in Fig.JIJd), MnPc adsorbs on 
Bi(110) with its Mn ion located directly on top of a dan- 
gling bond of the surface. This configuration suggests 
the formation of a local bond between MnPc and the 
dangling bond of the bismuth atom. 

The MnPc covered sample was subsequently annealed 
up to 130 K and exposed to CO gas partial pressure 
of 10 _T mbar for 90 s. After CO exposure the struc- 
ture of the molecular islands remains unchanged, but 
several molecules (~10%) exhibit now larger apparent 
height (Fig.HJb)). We identify these new molecules as 
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FIG. 1: Adsorption of MnPc and CO-MnPc on Bi(llO). 
(a) STM image of highly ordered MnPc island (7=0.2 nA, 
V=250mV, [44]). (b) MnPc island after exposure to CO 
(7=0.1 nA, y=180mV). CO-coordinated molecules can be 
distinguished by different apparent height. Inset images in 
(a) and (b) show schematic pictures of the chemical struc- 
ture of MnPc and CO-MnPc. (c) High resolution STM image 
of bare and CO-ligated MnPc (1=0.1 nA, V=-250mV). (d) 
Adsorption model of MnPc on Bi(110). 
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FIG. 2: (a) dl/dV spectra of MnPc and CO-coordinated 
MnPc in the bias range close to Ef, showing the zero bias 
anomaly. The dotted plots are the spectra measured on bare 
Bi surface as a reference, (b) Zero-bias feature of MnPc mea- 
sured at reduced temperature T=2.8K. The feature exhibits 
clear Fano line shape. Continuous line (red) corresponds to 
the fit by Fano equation. 

CO-ligated MnPc. Figure[TJc) shows a high resolution 
STM image of the new species embedded in the MnPc 
island. Similarly to MnPc, the CO-coordinated molecules 
exhibit fourfold symmetry, with a central protrusion ap- 
pearing 0.8 A higher than in case of MnPc, indicating 
that a single CO molecule bonds directly to the tran- 
sition metal ion. We find only singly coordinated CO- 
MnPc, contrary to previous studies reporting doubly co- 
ordinated porphyrin molecules jT3l [14] . 

The differential conductance (dl/dV) spectra mea- 
sured close to Ep reveals a pronounced anomaly at zero 
bias with a dip-like line shape on MnPc (Fig.[2ja)). 
We interpret it as a fingerprint of the Kondo effect, 
as it has been observed before for MnPc and other 



metal-phthalocyanines on surfaces [15-18 . The zero- 
bias anomaly exhibits a clear Fano line shape (Fig.[2ja) 
and (b)), commonly observed for Kondo ground states 
of magnetic adsorbates [19]. Other possible effect which 
could also cause a zero-bias anomaly, as e.g. inelastic spin 
flip excitation, have been carefully excluded by studying 
the temperature dependence of the line shape. A Kondo 
temperature of T^ nPc =22(±6) K can be approximated 
from a fit to the anomaly's line shape using the Fano 
equation f20ll2T]. 

When similar spectra is measured on CO-MnPc 
molecules, the ZBA appears much broader, with a 
slightly different lineshape (Fig. [2^ a)). The broadening 
of the zero-bias resonance indicates that the coordination 
to CO modifies the magnetic state of the molecule. This 
modification is reversible and can be controlled by selec- 
tive removal of CO molecules using STM tip. In order to 
detach an individual CO from its site, the STM tip was 
placed over CO-MnPc molecule, the feedback loop was 
opened and the sample bias was ramped. The resulting 
current- voltage characteristics is shown in Fig. [3^ a). A 
sudden drop of the tunneling current at a certain thresh- 
old voltage indicates detachment of the CO molecule 
from its coordination site (in most cases to a new site 
in a neighbor MnPc molecule). The resulting species re- 
covers the usual clove-like shape and a narrower zero- 
bias anomaly characteristic of a bare MnPc molecule 
(Fig.^c)). To establish the underlying mechanism of 
CO detachment we note that the threshold voltage in- 
creases linearly with the distance between the tip and 
the molecule, as shown in Fig.[3jb). This behavior indi- 
cates that the desorption of CO is induced by the electric 
field at the tunnel junction. A critical value of electric 
field e = l V/nm can be extracted from the slope of the 
line. 

The origin of the zero-bias anomaly on the CO-MnPc 
molecules cannot be deduced directly from our exper- 
iments. One could associate it to a fingerprint of a 
Kondo ground state with a larger Kondo temperature 
^co-MnPc about 50 (±io)K can be fitted). This is, 
however, unprovable considering the tt- acceptor charac- 
ter of the CO ligand. As a strong field ligand, CO affects 
the magnetic state of MnPc in two ways: (i) it increases 
the splitting of the <i-orbitals, and thus favors the low 
spin configuration in the complex [8 ; (ii) it withdraws 
the electron density from the bond to the ligand in the 
opposite coordination site (here the Bi dangling bond) 
and hence reduces the coupling of the d-orbitals to the 
substrate [2HI]. The latter effect, known in coordination 
chemistry as £rans-effect, should reduce the Kondo tem- 
perature, contrary to our experimental observation. 

To rationalize the experimental results we performed 
density functional theory (DFT) calculations of free and 
adsorbed MnPc and CO-MnPc molecules. All calcu- 
lations were performed using the Green's function for- 
malism implemented in our code ANT.G, which inter- 
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FIG. 3: Desorption of CO molecules from MnPc islands, 
(a) I-V characteristics recorded during desorption of CO. A 
sudden drop in current marks the threshold voltage for CO 
detachment, (b) Threshold voltage plot versus tip-sample dis- 
tance. Zero position of the tip corresponds to the feedback 
parameters V/=100mV and i/=100pA. Each point is an av- 
erage of over 30 measurements. Linear fit of the data indicates 
that the process is driven by electric field of 1 V/nm. (c) STM 
image and dl/dV spectrum of the molecule before and after 
controlled CO desorption. 

faces with Gaussian03/09 [22j [23] . We used a gen- 
eral gradient approximation (GGA) with the Becke and 
Perdew-Becke-Ernzerhof exchange-correlation functional 
[24] in combination with a common double local basis 
set LANL2DZ [25j[26 . Diff erent views of the computed 
system are shown in Fig.Qa). The Bi(110) surface was 
described using two bilayers embedded on a tight-binding 
Bethe lattice. Following the experimental results (Fig.[T]), 
MnPc and CO-MnPc molecules were placed with the 
manganese atom directly on top of a Bi atom with an op- 
timized distance to the surface and rotated seven degrees 
with respect to the high symmetry axis of the Bi(110) 
surface. In this configuration, four carbon atoms of the 
phthalocyanine cycle and the central Mn ion are in direct 
contact with Bi(110) dangling bonds. 

We simulate first the electronic and spin configura- 
tion of the free MnPc molecule, obtaining results in close 
agreement with previous studies 27 29 . The molecule 
has a spin state 5 = 3/2 due to three unpaired elec- 
trons localized on the manganese orbitals d xz /d yzi d z i 
and d xy . The d xz /d yz orbitals, hosting three electrons, 
are strongly hybridized with the phthalocyanine ring of 
the molecule, masking the atomic character of its local- 
ized spin. The d xy orbital hosts the most localized un- 
paired electron. 

On the Bi(110) surface, MnPc keeps two unpaired elec- 
trons in the d-orbitals. The net spin density in the d xz 
and d yz orbitals is reduced and results now from four 
occupied molecular orbitals with different atomic contri- 
bution depending on the spin orientation (Figj4^b)). In 
addition, a spin polarization arises in the ligands, op- 



posite to that in the Mn core. Thus, the total spin is 
reduced to S = 1, with two identifiable and localized un- 
paired electrons in the the d z i and d xy orbitals (Fig|4|[b)). 
We also find that a spin polarized density arises on Bi 
atoms directly underneath the adsorbed molecule, anti- 
ferromagnetically coupled to the spin on the Mn core and 
thus supporting the emergence of Kondo screening. 

The strength of a Kondo screening channel depends, 
to a first approximation, on the hybridization with the 
substrate. For MnPc on Bi(110), the spin-polarized local 
density (LDOS) of states projected on the Mn-d orbitals 
shows that each orbital hosting an unpair electron inter- 
acts differently with the surface (FigQd)). The Mn-d^ 
orbital hybridizes with a dangling bond of Bi and broad- 
ens substantially. The d xy orbital appears as narrow 
peaks in the PDOS spectrum. This indicates that this 
state remains highly decoupled from the surface, prob- 
ably to the absence of z-component. The result is that 
the two unpaired spins are screened via Kondo screen- 
ing channels with very different energy scale and, thus, 
Kondo temperature [17]. Such scheme suggests an un- 
derscreened spin if the experimental temperature lies be- 
tween the Kondo temperatures of each screening channels 

|23 ED- 

In the DFT framework, the orbital-substrate coupling 
is represented by a hybridization function which allows 
for a qualitative analysis of the feasible Kondo channels 
[32j[33]. The hybridization function is given by: X)d(^) = 
(cj + /i)I — Hjp — (Gd(^)) -1 , where uj and /i are the 
energy and the chemical potential respectively, I is the 
identity matrix, H^ s is the Kohn-Sham Hamiltonian and 
Gd is the one electron green's function projected onto d- 
subshell. Only the imaginary part of the hybridization 
function is relevant in this work (A = -7mSd(w)/7r). 
Specifically, the prerequisite for the Kondo effect of the 
host orbital having a finite coupling to the electron bath 
is equivalent to having non-zero value of the hybridiza- 
tion function A at the Fermi level [34 . The value of A 
calculated for different Mn d-orbitals of adsorbed MnPc 
is shown in Fig.Qg). The d z i orbital shows a peak in 
the hybridization function at the Fermi level indicating 
that this orbital is strongly screened by the substrate. On 
the contrary, the hybridization function of the d^-orbital 
has a zero value in a broad energy window in the vicinity 
of the Fermi level. These results suggest that the Fano 
resonance in the spectra of the MnPc molecule is a finger- 
print of a Kondo channel opened to screen the d z i orbital. 
The unpaired spin of the d^-orbital remains decoupled 
from the electron bath and this Kondo channel is closed. 
On the basis of these simulations, we conclude that the 
MnPc molecule retains part of its magnetic moment on 
the surface, lying in an underscreened many-body ground 
state. This is a consequence of the huge disparity in hy- 
bridization functions between spin channels, making this 
system unique in this regard. 

Next, we perform a similar analysis for the CO-MnPc 
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molecular complex. The CO molecule stays perpendic- 
ular to the MnPc molecular plane (Fig.Qa)). It is an- 
chored to the manganese ion via its carbon atom through 
a synergic 7r*-back-bonding, which strongly modifies the 
electron population and the spin of the Mn d-orbitals 
in the following manner. The d z i and the d xz /d yz or- 
bit als loose their magnetic moment: the former is emp- 
tied due to the overlap with a nonbonding orbital of CO 
and the latter is fully occupied due to the bond with the 
p x and the p y orbitals of CO. The spin-resolved LDOS of 
the adsorbed CO-MnPc projected on the Mn d-orbitals 
(FigQf)) confirms that the spin-polarization of the d z 2 
and d xz /d yz orbitals is reduced. Only the d^-orbital ap- 
pears to be relevant for the magnetism of the complex. 
The total spin of the MnPc-CO complex is reduced to 
S = 1/2, steaming from an unpaired electron localized in 
a molecular orbital with strong atomic d^-orbital char- 
acter but a finite contribution from the phthalocyanine 
cycle (see FigQc)). 

It thus appears clear that the modification of the ex- 
perimentally observed Kondo resonance upon CO coordi- 
nation is a direct consequence of the redistribution of the 
^-electrons of the Mn atom. However, the ZBA broad- 
ening cannot be understood as an increase in the Kondo 
screening because the responsible level is very localized. 
As seen in Fig.Qg), the d xy hybridization function is zero 
in the vicinity of the Fermi level, what confirms that also 
now no Kondo effect can be mediated by this orbital. 

This apparent inconsistency with the experimental ob- 
servations can be bypassed by noting that the PDOS on 
the d xy -orbital shows a highly localized branch pinned 
at the Fermi level. This state thus belongs to the low- 
est unoccupied orbital of the molecule (LUMO), and is 
the only responsible for the chemical potential line-up 
between the molecule and the substrate [35]. In these 
circumstances, charge fluctuations are expected to oc- 
cur at the d^y-orbital. When electronic correlations are 
considered, a mixed- valence regime is likely to emerge, 
as it is known for the /-electron compounds [34 . In 
this situation, one expects a resonance to appear in the 
differential conductance spectra at zero bias, similar to 
that in the Kondo regime but characterized by a larger 
width [36] [37]. It is interesting to note that the small 
hybridization with the surface prevents the finite spin in 
the d^-orbital from quenching, in contrast with the sce- 
nario reported in related systems [38] [39]. This allows 
for a many-body correlated state to emerge. 

Finally, we comment on the possible role of the sur- 
face states in the screening process on Bi(llO). As a 
semimetal, bismuth is characterized by a low density 
of bulk states at the Fermi level; its surfaces are how- 
ever strongly metallic due to the presence of the surface 
states. Therefore, the Kondo effect on bismuth can be 
expected to involve surface states rather than bulk states. 
However, strong spin-orbit coupling induces the spin po- 
larization of the surface states through the Rahsba ef- 
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FIG. 4: (a) Model of computed system: top view (left) and 
side view (right) of CO-MnPc (top) and MnPc (down). (b,c) 
Energy spectrum of adsorbed MnPc (b) and CO-MnPc (c), 
showing the relevant spin states. Bars in left (right) show 
the up(down) eigenvalues. The rate in parenthesis gives the 
atomic character of the molecular orbital, (d, f) Local density 
of states (PDOS) and (e,g) hybridization function (A) pro- 
jected on d-orbitals of adsorbed MnPc (d, e) and CO-MnPc 
(f, g). The dotted lines in (e) and (g) refer to the orbitals 
which are not relevant for the magnetism of the molecule. 



feet [TT]; the concept of screening in such system is not 
straightforward [40] . The Kondo effect involving the chi- 
ral spin-states has been recently intensively discussed in 
the context of topological insulators [4TH43] and can be 
understood assuming a complex structure of the Kondo 
cloud with a nontrivial spatial and spin dependence [43] . 

In summary, we have demonstrated that the magnetic 
ground state of individual metal-organic molecules can 
be reversibly modified by chemical coordination to ex- 
ternal molecular species acting as a ligands. Using a 
combination of low-energy scanning tunneling spectra 
and ab initio calculations, we find that the bonding of 
a CO molecule to the Mn 2+ center of MnPc molecules 
on Bi(110) leads to a substantial modification of its spin 
state. Our results indicate that the chemically induced 
magnetic transition corresponds to a change from a high- 
spin underscreened Kondo state to a mixed- valence state, 
where charge fluctuations destroy the spin 1/2 state of 
the CO-MnPc complex. Our results thus emphasize the 
ultimate limits of manipulating magnetic states of in- 
dividual molecules by using a combination of tunneling 
electrons and chemical control over the molecular spin, 
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which could facilitate molecular scale patterning of mag- 
netic motifs in molecular thin films. 
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